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Abstract 

We calculate the branching ratios and CP-violating asymmetries for B° — ► ton, con' , 4>n and 
4>n' decays in the perturbative QCD (pQCD) factorization approach. The pQCD predictions 
for the CP-averaged branching ratios are Br(B° — » cor/) = (2.7^ 1 ' ) x 10~ 7 , Br(B° — > con') = 
(0.75±g:|^) x 10~ 7 , and Br(B° -> 4>n) = (6.3±f;|) x 10~ 9 , Br(B° -> tfyrf) = (7.3±|;|) x 10~ 9 
which are consistent with currently available experimental upper limits. The inclusion of the 
gluonic contribution can change the branching ratios of B — > co((f))rf decays by about 10%. The 
direct CP-violating asymmetries for B° — » cor] and ton' decays are generally large in size. 
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I. INTRODUCTION 



As is well-known, the two-body charmless B meson decays provide a good place for 
testing the standard model (SM) and searching for the new physics signals. Among various 
B — > M\M.2 ( here Mj refers to the light pseudo-scalar or vector mesons ) decay channels, 
the decays involving the 77 or rj meson in the final state have been studied extensively 
during the past decade because of the so-called Krf puzzle or other special features. 

In Ref. for example, many decay modes involving 77 W meson were studied in the 
QCD factorization (QCDF) approach [2]. In the pQCD factorization approach [3|, on 
the other hand, the B — > Krj('\ prj^'\ tct]^ and ry^ry^' decays have been calculated in 
Refs. 0, H, @, 0]. The B s —>■ (jr, p, 00, <p)i]^ decays, furthermore, have also been studied in 
the pQCD approach [gj]. 

In this paper, we will perform the leading order pQCD calculation for B — > urj^ and 
(f)7]^ decays. Besides the usual factorizable contributions, we here are able to evaluate 
the non-factorizable and the annihilation contributions as well. Besides the dominant 
contributions from the qq components of 77''' mesons, the contribution from the possible 
gluonic component of rf'^ meson will also be included by using the formulae as given in 
Ref. @. 

On the experimental side, only the upper limits on the branching ratios of B — > u((f>)7)('' 



decays are available now [10, 111 



Br(B° -> (f)T)) < 0.6 x 10~ 6 , Br(B° -> <frrf) < 1.0 x 10~ 6 . (1) 
Br(B° -»■ lot]) < 1.9 x 10~ 6 , Br(B° -> urf) < 2.8 x 10~ 6 . (2) 

But these decays could be measured with good precision in the forthcoming LHC-b ex- 
periments if their decay rates are larger than 10 -8 . 

For B — > u((j))r)^ decays considered here, it is generally believed that the qq component 
of r)^ meson provide the dominant contribution, but it is still very difficult to calculate 
reliably the gluonic contribution from the possible glunoic component of 77 W meson 
Of course, great efforts have been made for this problem and some progress have been 
achieved recently 0, [l2| in order to explain the large B — > ^77' decay rates. In Ref. B , 
for instance, the authors found that the possible gluonic contribution is small for both 
B — > 77 and B — > rj' form factors 0, 13| . 

This paper is organized as follows. In Sec. [Til we gi ye a brief review for the PQCD fac- 
torization approach. In Sec. IHH we calculate analytically the related Feynman diagrams 
and present the various decay amplitudes for the studied decay modes. In Sec. [TVl we 
show the numerical results for the branching ratios and CP asymmetries of B — > u{<f))rf'' 
decays. The summary and some discussions are included in the final section. 



II. THE THEORETICAL FRAMEWORK 

At present, there exist two popular factorization approaches to calculate the hadronic 
matrix element < MiM 2 \Oi\B >: the QCDF approach [2| and the pQCD approach 
(3, 14, lE)]. The pQCD approach has been developed earlier from the QCD hard-scattering 
approach jl5|. Some elements of this approach are also present in the QCD factorization 
approach The two major differences between these two approaches are (a) the 
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form factors are calculable perturbatively in pQCD approach, but taken as the input 
parameters extracted from other experimental measurements in the QCDF approach; and 
(b) the annihilation contributions are calculable and play an important role in producing 
CP violation for the considered decay modes in pQCD approach, but it could not be 
evaluated reliably in QCDF approach. Of course, the assumptions behind the pQCD 
approach, specifically the possibility to calculate the form factors perturbatively, are still 
under discussion [l6l |. More efforts are needed to clarify these problems. 

In pQCD approach, the decay amplitude is separated into soft ($), hard(H), and 
harder (C) dynamics characterized by different energy scales (t, m b , M w ). It is conceptu- 
ally written as the convolution, 

A(B -> M X M 2 ) ~ Jd%d%d% Tr [C(t)^ B (h)^ Ml (k2)^M 2 {h)H{h, k 2 ,k 3 ,t)\ , (3) 

where the hi are the momenta of the light quarks included in each of the mesons, and 
Tr denotes the trace over Dirac and color indices. C(t) is the Wilson coefficient which 
results from the radiative corrections at short distance. In the above convolution, C(t) 
includes the harder dynamics at larger scale than Mb scale and describes the evolution 
of local 4-Fermi operators from my/ (the W boson mass) down to t ~ scale, 
where A = M B — vrib- The function H(k\, k 2 , k 3 , t) is the hard part and can be calculated 
perturbatively. The function <3>m is the wave function which describes hadronization of 
the quark and anti-quark to the meson M. While the function H depends on the process 
considered, the wave function <E> M is independent of the specific decay process. Using the 
wave functions determined from other well measured processes, one can make quantitative 
predictions here. 

Since the b quark is rather heavy, we consider the B meson at rest for simplicity. It is 
convenient to use light-cone coordinate (p + ,p~,Pr) to describe the meson's momenta, 

P ± = -^(P°±P 3 ), and p T = (p\p 2 ). (4) 

Using these coordinates the B meson and the two final state meson momenta can be 
written as 



Mr Mr o Mr 



respectively, where = m ^^/ 'm B ; and the terms proportional to the ratio m 2 ^ {l) /m B 
have been neglected. 

For the B un decay considered here, only the u meson's longitudinal part con- 
tributes to the decay, its polar vector is = J^^ (1, — r^, T )- Putting the light (anti-) 
quark momenta in the B, uj and rj meson k\, k 2 , and k^, respectively, we can choose 

h = (x 1 P+,0,k 1 T), k 2 = (x 2 P 2 + ,0,k 2T ), k 3 = (0,x 3 P 3 -,k3 T ). (6) 

Then, the integration over k~[ , k 2 , and k% in Eq.Q will lead to 



A(B — ► Lurj) ~ / dxidx 2 dx s bidbib 2 db 2 b 3 db 3 

•Tr [Cm B (x 1 ,b 1 )^(x 2 ,b 2 )%(x 3 ,b 3 )H(x l ,b l ,t)S t (x t )e- s ^] , (7) 
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where 6j is the conjugate space coordinate of kir, and t is the largest energy scale in the 
function H{x^ 6j, t). The large logarithms ln(m^/t) are included in the Wilson coefficients 
C(t). The large double logarithms (In 2 Xj) on the longitudinal direction are summed by the 
threshold resummation, and they lead to the function St(xi) which smears the end-point 
singularities on X{. The last term, e -5 ^-*, is the Sudakov form factor which suppresses the 
soft dynamics effectively 1?J • Thus it makes the perturbative calculation of the hard part 
H applicable at an intermediate scale, i.e., the Mb scale. We will calculate analytically 
the function H(xi, b i: t) for the considered decays in the first order in an a s expansion and 
give the convoluted amplitudes in the next section. 



Wilson Coefficients 



For the B° -> 
can be written as 
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B° — > (f)T]('' decays, the related weak effective Hamiltonian H, 



eff 



U 



eff 



G f 

V2 
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v ub v: d {cmo u m + cmo u m) - v tb v; d £ c^) o^) 



i=3 



where Cj(/i) are Wilson coefficients at the renormalization scale /i and Oi are the four- 
fermion operators. For b — > d transition, for example, these operators can be written 

as 
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d_ a j^Lu a ■ Ufs'J^Lbp , 
da^Lbp ■ J2 q > q'pl^Lq'a , 

\d a ^Lb p ■ Y^ q , e ql q'^^Rq' a 

\d a YLb p ■ J2 q > a q 'q'piM'a 



(9) 



where a and j3 are the SU(3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), R = (1 + 75). The sum over q' runs over the 
quark fields that are active at the scale /z = 0(m b ), i.e., q'e{u, d, s, c, b}. For the Wilson 
coefficients Cj(/x) (i — 1, . . . , 10), we will use the leading order (LO) expressions, although 



the next-to-leading order (NLO) results already exist in the literature [18|. This is 
the consistent way to cancel the explicit fi dependence in the theoretical formulae. For 
the renormalization group evolution of the Wilson coefficients from higher scale to lower 
scale, we use the formulae as given in Ref. directly. At the high m w scale, the leading 
order Wilson coefficients Ci(M\y) are simple and can be found easily in Ref. 18]. In 
the pQCD approach, the scale t may be larger or smaller than the mj scale. For the 
case of nib < t < mw, we evaluate the Wilson coefficients at t scale using the leading 



logarithm running equations, as given in Eq.(Cl) of Ref. jl9[. For the case of t < nib, we 
then evaluate the Wilson coefficients at t scale by using Cj(m&) as input and the formulae 
given in Appendix D of Ref. 1^ • 
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B. <j> — u> mixing and rj — rf mixing 



For the physical isoscalars cj) and ui meson, we use the "ideal" mixing scheme [20 



0(1020) = -ss, 



1 
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—j= [uu + dd\ 
v2 



(10) 



Such ideal mixing are supported by the known experimental measurements [20 . 

For the rj — rj' system, there exist two popular mixing basis: the octet-singlet basis and 
the quark-flavor basis [2l|, |22| . Here we use the quark-flavor basis 2l| and define 

t] q = (uu + dd)j\p2, T] s = ss. (11) 
The physical states rj and rf are related to rj q and r) s through a single mixing angle <fi, 



U{<j>) 



Vs 



cos . 
sin ( 



— sin i 

cos <b 



Vs 



The corresponding decay constants f q , f s , f^' s and f^f have been defined in Ref. [21 



(12) 



as 



< 0\s^J 5 s\r] s {P) > 



V2 Jq ' 



<0\qYl,q\v ( '\P)> = -^=f 9 M P" 



< 0\sYl5s\v [ '\P) > 



V2 



(13) 



(14) 



while the decay constants f%' s and fi s are related to /„ and f s via the same mixing matrix, 



fq fs 

J r\ J r\ 

fl fs 

Jn' Jn> 



u{4>) 



o fs 



(15) 



The three input parameters f q , /«, and in the quark-flavor basis have been extracted 
from various related experiments 21, 22| 



f q = (1.07 ± 0.02)/*, f s = (1.34 ± 0.06)/*, 



39.3° ± 1.0°, 



(16) 



where f n = 130 MeV. In the numerical calculations, we will use these mixing parameters 
as inputs. 

Although rj and r/ are generally considered as a linear combination of light quark pairs 
uu, dd and ss, it should be noted that a gluonic content of rf meson may bejieeded to 

231. For 



interpret the anomalously large branching ratios of B — > Kr/' and J/^f 
B — > (p, 7r)ryW decays, however, the good agreement between the pQCD predictions @] 
and the measured values for their branching ratios suggest that the gluonic contribution 
of 7]^ meson should be small. In Ref. 0], very recently, the authors calculated the flavor- 
singlet contribution to the B — > 77 ^ transition form factors from the gluonic content of 
the 7/') meson. They found that the gluonic contribution is small for both B — > 77 and 
B — > rf form factors. 
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Of course, more studies are needed to provide a reliable and precision calculation for 
the gluonic contribution to the final state involving 77 W meson. For the sake of simplicity, 
we here firstly neglect the possible gluonic content of 77 and rj meson, and use the quark- 
flavor mixing scheme for 77 and 77' meson as defined in Eq. (TP2"1) . We will calculate the 
effects of a non-zero gluonic admixture of 77' in next section, and treat them as one kind 
of theoretical uncertainties. 



C. Wave Functions 

Now we present the wave functions to be used in the integration. For the wave function 
of the heavy B meson, we take 



B 



V2N C 



i B + M B )7 5 jB (k 1 



(17) 



Here only the contribution of the Lorentz structure 0s(ki) is taken into account, since 
the contribution of the second Lorentz structure 0b is numerically small 0, 25] and has 
been neglected. For the distribution amplitude 0b in Eq- fTTTj) . we adopt the model 



4>b(%, b) = Nbx 2 (1 — x) 2 exp 



M 2 x 2 



2u 



2 -TjMY 



where uj^ is a free parameter and we take u>b = 0.4 ± 0.04 GeV in numerical calculations, 
and Nb = 91.745 is the normalization factor for u>b = 0.4. This is the same wave functions 
as being used in Refs. [19|, |24| , which is a best fit for most of the measured hadronic B 
decays. 

The wave function for dd components of 77 ^ meson is given by 



_H5_ 



Vdd 



X 



Vdd 



x 



+ (m 



ids 




(19) 



where p and x are the momentum and the momentum fraction of r] d ^ respectively, while 
0^ _, 0^ _ and _ represent the axial vector, pseudoscalar and tensor components of the 

'dd 'dd 'dd I — 1 

wave function respectively. Following Ref. [4|, we here also assume that the wave function 
of 77 d( j is same as the tt wave function based on SU(3) flavor symmetry. The parameter ( 
is either +1 or —1 depending on the assignment of the momentum fraction x. We also 
assume that the wave function of the uu component is the same as the wave function of 
dd based on the isospin symmetry between the up and down quark. For the wave function 
of the ss component, we also use the same form as given in Eq. (fT9l for dd component 
but with different chiral enhancement and some other changes to be specified later. 

The explicit expressions of chiral enhancement scales m, 
are given by (§] 



m. 



Vdd 



m 



and tUq = m, 



m 



gq 



2m n 



2m , 



r 2 2 1 1 2-2 
- — cos + m, sin - 

2m n 1 1 



r 2 2 j 1 2 • 2 
■ [ml, cos <p + m v sin 



2m. 



u 
u 

V2f s 



(m 2 , — m 2 ) cos0si?7^ 



(m 2 , — m 2 ) cos 0sm0], 



(20) 
(21) 
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and numerically 



m, 



o 



1.07MeV, 



1.92GeV 



(22) 



for m v = 547.5 MeV, Try = 957.8 MeV, f q = 1.07/,,, f s = 1.34/ w and = 39.3°. 

For the distribution amplitude 0^ , 0f and 0^ , we utilize the results for tt meson 
obtained from the light-cone sum rule 26J including twist-3 contributions: 

.A , . 3 



Vg(s) 



[X 



y [21(1 -2x) 4 - 14(1 -2x) 2 + l]j, 



(23) 



0^ (x) 



1 



2./2N~ c JHK °' [ 2 

1 / 27 2 81 2 

-3%^3 - ™P n „^ - 



2 rj 7 9 ( 3 ) 



) [3(1 - 2xf - 1] 



Vq(s) 



[X 



[35(1 - 2x) 4 - 30(1 - 2xf + 3] } , 

o 

=/,(,)(! - 2x) 



1 



+ (5r? 3 - -%^3 - — p 2 , --pl a^HlO^-lOz + l) 
6 2 20 Vq{s) 5 



(24) 



(25) 



with the updated Gegenbauer moments [27], |28 



"2 



0.115, a 



V q (s) 



-0.015, p v = 2m q /m. 



qqi 



p Vs = 2m s /m ss , r/3 = 0.015, uj 3 = -3.0. 



(26) 



For 5 — > wryW an( } 0^(0 decays, only the longitudinal polarized component $^ for 
V = (cu, 0) will contribute. The wave function for example, can be written as 



[m w /l0 w (x)+ / l #^(x) +m w /0^(x)] 



(27) 



where the and p is the polarization vector and the momentum of the u meson, the first 
term in above equation is the leading twist wave function (twist-2), while the second and 
third terms are subleading twist (twist-3) wave functions. For the case of V — 0, its wave 
function is the same in structure as that defined in Eq. (127]) . but with different mass and 
distribution amplitudes. 

For the light meson wave function, we neglect the b dependent part, which is not 
important in numerical analysis. The distribution amplitudes <f><j>(x) and 0^ s (x) are given 
In 1} 



V6 f 

J<t> 

2~7I 

3 

WE 



x(l — x), 

3(1 - 2x) 2 + 1.68C 4 1/2 (1 - 2x) + 0.69 



l + (l-2x) In 



x 



X 



J(t 



3(1 - 2a;) (4.5 - ll.2x + 11.2x 2 ) + 1.38 In 



x 



X 



(28) 
(29) 
(30) 
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For the u meson wave function, we use 30, 31 



<j> u (x) = ^=f u x{l -x)\l + 0.18O 2 3/2 (l - 2x)l , (31) 

= J*L {3(1 - 2xf + 0.3(1 - 2x) 2 [5(1 - 2xf - 3] 

+0.21[3 - 30(1 - 2x) 2 + 35(1 - 2x) 4 ]} , (32) 

C(x) = - 2x) [1 + 0.76(10x 2 - 10a: + 1)] . (33) 

The relevant Gegenbauer polynomials are defined by [29| 

C y\t) = \ (5t 2 - 1) , (34) 

C\ /2 (t) = I (35t 4 - 30t 2 + 3) . (35) 



III. PERTURB ATIVE CALCULATIONS 

In this section, we will calculate and show the decay amplitude for each diagram 
including wave functions. The hard part H(t) involves the four quark operators and the 
necessary hard gluon connecting the four quark operator and the spectator quark. We 
first consider B — > ujt] decay, and then extend the calculation to other decay modes. 
Analogous to the B — > prj® decays in (^], there are also eight type diagrams contributing 
to B — > wq and ujrf decay, as illustrated in Figure 1. 

For B — > wq decay, we firstly consider the usual factorizable diagrams 1(a) and 1(b). 
The operators Oi,2 and 03,4,9,10 are (V — A)(V — A) currents; the sum of their amplitudes 
is given by 



F P = 



4v / 27rGi?(7 F m 4 3 / dx\dx 3 \ 61^6163^63 <f>s(xi, 61) 
Jo Jo 

X { [(1 + X 3 )<^(x 3 , 63) + (1 - 2X3)^(0^X3, 63) + f v 63))] 

■a s {t\) h e (x lt x 3 , 61, 63) exp[-Sab(tl)) 
+2r v 4>^(x 3 , b 3 )a s (t 2 e )h e (x 3 , xi, 63, 61) exp[-S ab (t 2 e )]} . (36) 

where the ratio r v = m^/ % with m 1 ^ = ml or m s as defined in Eqs. (120112 ip ; the color 
factor Cf = 4/3, 4>b and 4>^' P ' T are the light-cone distribution amplitudes (LCDAs) of the 
heavy B meson and the light 77 meson. The functions h e , the scales t\ and the Sudakov 
factors S a b will be given explicitly in Appendix IA1 

The operators 05,6,7,8 have the structure of (V — A)(V + A). In some decay channels, 
some of these operators contribute to the decay amplitude in a factorizable way. Since 
only the vector part of the (V + A) current contributes to the vector meson production, 
(r}\V-A\B)(u\V + A\0) = (ij\V- A\B)(u\V- A\0), that is 



1P1 



(37) 
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(g) (h) 

FIG. 1: Feynman diagrams contributing to the B — * ujrf 1 ^ decays when the rj^ meson picks up 
the spectator quark, where diagram (a) and (b) contribute to the B — > rf 1 ^ form factor . 

For other cases, one needs to do a Fierz transformation for the corresponding operators to 
get the right flavor and color structure for factorization to work. We may get (S+P) (S—P) 
operators from (V — A){V + A) ones. Because neither scalar nor the pseudo-scalar density 
give contribution to a vector meson production, (u\S + P\0) = 0, we get F e P2 = 0. 

For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 63 can be performed using the S function 5(6 3 — b\), leaving 
only integration of 61 and 62- The decay amplitude for (V — A)(V — A) and (V — A)(V+A) 
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operators can be written as 

M e = -M[ 2 = —=G F irC F m B / dx\dx 2 dx 3 / bidbib 2 db 2 0s(xi, fa)(f) u; (x 2 , fa) 
v3 Jo Jo 

x { [2x 3 r v (f)T(x 3 , fa) - x 3 (p^(x 3 , fa)] 

■a s (t f )h f (xi, x 2 , x 3 , fa, fa) exp[-S cd (t f )]}, (38) 

32 f 1 f°° 

Mf 1 = -=GF^CF'Tu)'m A B I dx\dx 2 dx 3 \ fadfab 2 db 2 (f)B{xi,bi) 

v3 Jo Jo 

■ { [x 2 (j)^(x 3 , fa) ((f) s U! (x 2 , fa) - <pl(x 2 , fa)) + r v ((x 2 + x 3 ) ((^{x 3 , fa) 
•Glfa, fa) + 4%(x 3 , fa)(j)l(x 2 , fa)) + (x 3 - x 2 ) (</>%(x 3 , fa)(f) t ijJ (x 2 , fa) 
+(j)T(x 3 ,fa)(j) s U} (x 2 ,fa)))] a s (t f )h f (x 1 ,x 2 ,x 3 ,fa,b 2 ) exp[-S cd (t f )]} . (39) 

For the non-factorizable annihilation diagrams 1(e) and 1(f), there are three kinds of 
decay amplitudes. For the (V — A)(V — A) operators, the decay amplitude M a is 

16 f 1 f°° 

M a = —=kG pC F m B / dx\dx 2 dx 3 I fadfab 2 dfa <p B ( x i,fa) 
v3 Jo Jo 

■ { [r u r v (x 3 - x 2 ) [tf(x 3 , fa)4>l{x 2 , fa) + 0^(0:3, fa) 0* (x 2 , fa)] 

+W v (x 2 + x 3 ) [<fr^(x 3 , b 2 )(j) s 0J (x 2 , fa) + 0j(x 3 , fa)(j)l(x 2 , fa)] 

+x 3 (p UJ (x 2 , fa)(p^(x 3 , fa)] ■ a s {t))tif(xi, x 2 , x 3 , fa, fa) exp[-S ef (t 4 f )] 

- [r u r v (x 2 - x 3 ) [4>^{x 3 , b 2 )^{x 2 , fa) + <p^(x 3 , fa)<p s w (x 2 , fa)] 

+r UJ r v [(2 + x 2 + x 3 )<p^(x 3 , fa)(f)l(x 2 , fa) - (2 - x 2 - x 3 )(f^(x 3 , b 2 )^{x 2 , fa)] 

+x 2 (p UJ (x 2 ,fa)(f)^(x 3 ,fa)] ■a s (t 3 f )h 3 f {x 1 ,x 2 ,x 3 ,fa,fa)exp[-S ef (t 3 f )]} . (40) 

For the (V - A)(V + A) and (S - P)(S + P) operators, we have 

16 /*°° 
M[ l = —=GF^CFfnB \ dx\dx 2 dx 3 I fadfafadfa4>B{xi,fa) 
v3 Jo Jo 

■ { [^2^0^x3, fa) {4> s w {x 2 , fa) + (j)i(x 2 , fa)) - x 3 r v ((f)^(x 3 , fa) + (j%(x 3 , fa)) 
•<f>u,( x 2, fa)] ■ a s {t))h){xi, x 2 , x 3 , fa, fa) exp[-S ef (t 4 )] 

+ [(2 - x 2 )r^(x 3 , fa) (<f>' u {x 2 , fa) + (f)l(x 2 , fa)) - (2 - x 3 )r v (tf(x 3 , fa) 

+(j) T n {x 3 ,fa)) <pu{x 2 ,fa)] ■ a s (t 3 )h 3 f (x 1 ,x 2 ,x 3 ,fa,fa)exp[-S ef (t 3 f )]} , (41) 

16 f 1 f°° 

M^ 2 = —TrGpCFm^ / dxidx 2 dx 3 / fadfab 2 dfa(f) B (xi,fa) 

v3 Jo Jo 

■ { [x 2 ^{x 2 , fa)^(x 3 , fa) + r u r v ((x 2 - x 3 ) (0j(x 3 , fa)(f>l J (x 2 , fa) + <p T q (x 3 , fa) 

■<PI{X 2 , fa)) + (X 2 + X 3 ) (0^X3, fa)^(x 2 , fa) + <fi(x 3 , 6 2 )0L(^2, fa)))] 

•a s (t 4 f )h 4 f (x 1 ,x 2 , x 3 , fa, fa) exp[-S ef (t 4 f )] 

- [xz<pu,{x 2 , fa)<p^(x 3 , fa) + r^r^ {(x 3 - x 2 ) (<f>^(x 3 , fa)4> t ul (x 2 , fa) + 0^(x 3 , fa) 
4l{x 2 , fa)) + (2 + x 2 + x 3 )tf(x 3 , fa)<j)-l(x 2 , fa) -(2-x 2 - x 3 )(^(x 3 , fa) 
■<t>l{x 2 ,b 2 ))] -a^hjix^x^x^b^b^expi-Sefit 3 )]} . (42) 

10 



The factorizable annihilation diagrams 1(g) and 1(h) involve only rj and u wave func- 
tions. There are also three kinds of decay amplitudes: F a is for (V — A)(V — A) type 
operators, F PI and F^ 2 is for (V — A)(V + A) and (S — P)(S + P) type operators, 
respectively: 

F = -F P1 

pi pOO 

= -A\/2G F -nC F m% I dx 2 dx 3 / b 2 db 2 b 3 db 3 { [x 3 (j)^(x 3 , 63)^(^2, b 2 ) 
Jo Jo 

+2r w r v ((x 3 + l)rf(x 3 , b 3 ) + (x 3 - 1)0* (x 3 , 6 3 ))C(x 2 , 62)] 

•a S (^e)^a(^2, ^3, &2, 63) exp [S gh (tj!)] 

- [a;20^(a;3,&3)0^(x2,& 2 ) 

+2ry^((x 2 + 1)0* (x 2 , 6 2 ) + (x 2 - 1)^(^2, b 2 ))^{x 3 , 63)] 

•a s (tg)/i a (x 3 ,x 2 ,6 3 ,6 2 )exp[-S' g / l (tg)]} , (43) 
/» 1 ^00 
F^ 2 = —%\plG p-nC pim^ j dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

x { [^(^fe, 63) - 0j(»3, h))^(x 2 , 62) + 2r w 0^(x 3 , 6 3 )0* (x 2 , b 2 )] 

■oi s {tl)h a (x 2 , x 3 , b 2 , 63) exp[-S' 9/l (^)] 
+ [^2^(0^,(^2, 62) - <Pl(x 2 , 6 2 ))0^(x 3 , 63) + 2r J? w (x 2 , b 2 )<p^(x 3 , b 3 )] 

■a s (tj)h a (x 3 ,x 2 ,b 3 ,b 2 )exp[-S gh (tf)]} . (44) 

In the above equations, we have assumed that x\ « (x 2 ,x 3 ). Since the light quark mo- 
mentum fraction x± in the B meson is peaked at the small region, while quark momentum 
fraction x 3 of rj is peaked around 0.5, this is not a bad approximation. The numerical 
results also show that this approximation makes very little difference in the final result. 

By exchanging the position of 77 W and u in Fig. (JTJ, one finds the new Feynman 
diagrams are illustrated in Fig. (T5]). Just like the case of Fig. (JTJ), we firstly consider the 
factorizable diagrams 2(a) and 2(b). The decay amplitude F e induced by inserting the 
(V — A){V — A) operators is 

pi poo 

F eLU = 4V2G -p7i -CfTtib / dx\dx 3 \ bidbib 3 db 3 <pB{xi,bi) 

Jo Jo 

• { [(1 + x 3 )0 w (x 3 , 63) + (1 - 2x 3 )r w (0* (x 3 , 63) + 0L(^3, 63))] 
■a s (tl)h e (xx, x 3 , 61, fo 3 ) exp[-S ab (tl)} 

+2r u) <p s u) {x 3n b 3 ) ■ a s (t 2 e )h e (x 3 ,x 1 ,b 3 ,b 1 )exp[-S ab (t 2 e )}} , (45) 

The Fig. 2(a) and 2(b) are also the diagrams being used to extract the form factor Aq^. 
According the definition in Ref. 25j, we can extract Aq w from Eq. (I45|) . 
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FIG. 2: Feynman diagrams contributing to the B — > wr/W decays when the w boson picks up 
the spectator quark. 



From Fig. 2(a) and 2(b), it is easy to find the decay amplitude F^} and F^j: 



(46) 



/>1 />oo 

F^u = %\f2GpKCpr ri m B \ dx\dx 3 / bidbib 3 db 3 (j)B{xi,bi) 

Jo Jo 

■ { [Ms* fo s) + r u {(x 3 + 2)C(^3, 63) - x 3 (Pt(x 3 , 63))] 
■a s {t\)h e {x x , x 3 , 61, 63) exp [-5^(^)1 

+ (xi0 w (x 3 , 63) + 2r w 0" (ar 3 , 63)) a s (^)/i e (^ 3 , 21, 63, 61) exjp[-S ab (t 2 e ))} . (47) 
For the non-factorizable diagrams 2(c) and 2(d), the corresponding decay amplitudes 



are 



16 /*°° 

M euJ = —G F ixC F m B / dx 1 dx 2 dx 3 / b 1 db 1 b 2 db 2 (f)B{xi,b 1 )(f)^(x2,b2) 

v3 Jo Jo 

• {^3 [<pu{x3,bi) ~ 2^0^x3,61)] a 3 (tf)hf(x 1 ,X2,x 3 ,b 1 ,b 2 )exp[-Scd{tf)]} ,(46 
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and 



AC = 0, M!s M eu) . (49) 

For the non-factorizable annihilation diagrams 2(e) and 2(f), again all three wave 
functions are involved, the three kinds of decay amplitudes are of the form 

16 f 1 f 00 

M auJ = —=Gf-kCfV(i b I dx\dx 2 dx 3 I b\dbib 2 db 2 4>b{x\, 61) 
v3 Jo Jo 

■ { [x 3 (t>oj(x 3 , 6 2 )0^(x 2 , 6 2 ) + ?vr T; ((x 3 - x 2 ) {(/>%(x2, 62)0^x3, 6 2 ) + 4\ j {x 2 , 6 2 ) 
■(f>t(x 3 , 6 2 )) + (x 3 + x 2 ) (rf{x 2 , b 2 )(f)l(x 3 , 6 2 ) + 0j(x 2 , b 2 )<f)l{x 3 , 6 2 )))] 
■a s (tj)h 4 f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S ef (t 4 f )] - [x 2 <f> u (x 3 , b 2 )<f)^(x 2 , b 2 ) 
+r u r v ((x 2 - x 3 ) (rf(x 2 , 62)0^X3, 6 2 ) + <fr^(x 2 , b 2 )(p s UJ (x 3 , 6 2 )) + r u r ri 

■ ((2 + x 2 + x 3 )tf(x 2 , 62)0^x3, 6 2 ) -(2-x 2 - x 3 )(j)T(x 2 , b 2 )(t>i(x 3 , 6 2 )))] 
■a s (t 3 f )h 3 f (x u x 2 ,x 3 ,b u b 2 )enp[-S ef (t 3 f )}} , (50) 

16 /*°° 

= -j^GF^Cpm^ I dx\dx 2 dx 3 \ b\db\b 2 db 2 0s(:ri,6 1 ) 

v3 Jo Jo 

■ { [xir^^x^ 6 2 ) (rf(x 2 , 6 2 ) + 4%(x 2 , 6 2 )) - x 3 r w (0* (x 3 , 6 2 ) + <f>l(x 3 , 6 2 )) 
•0^(^2, 6 2 )] • a s (t))h)(x u x 2 , x 3 , 61, 6 2 ) exp[-S ef (t 4 f )] 

+ [(2 - x 2 )r v <i> u (x3, 6 2 ) (<j)^(x 2 , b 2 ) + </)^(x 2 , 6 2 )) - (2 - x 3 )r u (0* (x 3 , 6 2 ) 
+0L(^3, & 2 )) 0^(x 2 , 6 2 )] • a s (tj)^(xi, x 2 , x 3 , 61, 6 2 ) exp[-S e/ (^)] } , (51) 

16 /*oo 

= — j=tiG pC I dx\dx 2 dx 3 \ bidbib 2 db 2 4>B{xi,bi) 
v3 Jo Jo 

• { [x 2 w (a;3, b 2 )(j)^(x 2 , 6 2 ) + r^r^ ((x 2 - x 3 ) (rf(x 2 , b 2 )(f)l(x 3 , b 2 ) + $(x 2 , b 2 ) 
•<#L(*3, 62)) + (x 2 + x 3 ) (^(x 2 , 62)0^x3, 6 2 ) + 0^(x 2 , 6 2 )0^(x 3 , 6 2 )))] 
•a s (t 4 f )tif(x 1 ,x 2 , x 3 , 61, 6 2 ) exp[-S ef (t 4 f )] 

- [x 3 (f)u 1 {x 3) b 2 )(p^(x 2 , 6 2 ) + ((x 3 - x 2 ) (0j(x2, 62)0^x3, 6 2 ) + 0^(0:2, 6 2 ) 
•0^3, 62)) + (2 + x 2 + x 3 )0j'(a;2, 6 2 )0* (x 3 , 6 2 ) - (2 - x 2 - x 3 )4%(x 2 , b 2 ) 
•^3,62))] •a s (^)^(x 1 ,x 2 ,x 3 ,6 1 ,6 2 )exp[- ( S e/ (^)]}. (52) 

For the factorizable annihilation diagrams 2(g) and 2(h), we have 

F a . = -F a , F£ = -F?, F£ = -F?. (53) 

Following the same procedure, it is straightforward to calculate the decay amplitudes 
for the B — > 0r/W decays. It is worth of mentioning that all the eight Feynman diagrams in 
Fig. 1 after the replacements of the u meson by the meson will contribute to B — > 0r/') 
decays; but Figs.2(a)-2(d) could not contribute to the same decay since = —ss in the 
ideal mixing scheme of u — system. 
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Combining the contributions from the different diagrams, the total decay amplitude 
for B° — > u>7] can be written as 

V2M(lu V ) = FeF^U 



in id + -C 2 

<t ( \c, + jjc 4 + 2C 5 + jjc 6 + l -C 7 + l -C s + \c, - l -C lG 



+M e F 1 (<f ) ) 
+F 



L,^ 2 - Li ( C 3 + 2C 4 — 2C 6 — -C 8 — -Cg + 



( Cx + -C 2 ) I 



~it (^3 + ^4 - 2C 5 - jjc 6 - ic 7 - ic 8 + ic 9 - ic 10 J /, d 

— Ct ( C*3 + ^C*4 — Cg — 7^C 6 + TJ^V + gC*8 — ^Cg — r^IO ) 



Zt[^C s + Ce-~C 7 -zC s )tf 



6 



+F P2 

+M eu) F 1 (4>) 
+M eu F 2 (<p) 
+ (M a + M auJ )F 1 (<f ) ) 



1 

-( 

6 



£«C 2 — ^ ( C*3 + 2C 4 + 2C 6 + -Cs — -Cg + -jCiQ 
— Li [ C \ -rCfi- -Cg — 2^ 10 



£«C 2 — it ( C 3 + 2C 4 — —Cg 



-C 



10 



- (M* + Mf 1 + M£ ) ^(0)6 (c 5 - ^C 7 ) 



(Mf 2 +M ( S)F 1 (0)e t (2C 6 + -C S 



where ^ = V* b V ud , it = V^V td , and 

FM 



COS 



F 2 {4>) = - sin. 



(54) 



(55) 



are the mixing factors. The Wilson coefficients Cj = Cj(t) in Eq. (I54"|) should be calculated 
at the appropriate scale t using equations as given in the Appendices of Ref. [19| . 



By doing the replacements of f u 
amplitude for B° — > <fir) decay: 



and 



kA,PT 



L A,P,T 



one obtains the decay 



M{<h) = F e i t \C 3 + l -C± + C 5 + l -C, - l -C 7 - l -C % - l -C 9 - ic 10 J F x (0) 

+M e i t (c 4 - C 6 + ^C 8 - ic 10 J F x (0) + (M a + M^) ^ (c 4 - ic 10 J F 2 (<f>) 
+ (Mf 2 + M^) i t (c 6 - l -C^j F 2 (0). (56) 
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The complete decay amplitude Ai(urj') and Ai(<f)f]') can be obtained easily from 
Eg. (|54|) and Eq. (|56p by the following replacements 

rd rs ed rs 

Jr\i J rj Jrj'i Jrj'i 



>/2' 

F 2 (0) — ^(0) = cos0. (57) 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will present the numerical results of the branching ratios and CP 
violating asymmetries for those considered decay modes. 



A. Input parameters 

We show here the input parameters to be used in the numerical calculations. 
The masses, decay constants, QCD scale and B° meson lifetime are 

A^T 4) = 250MeV, / w = 130MeV, f B = 190MeV, 
U = 237MeV, /J = 220MeV, M w = 80.41GeV, 
f u = 195MeV, fl = 140MeV, = 1.02GeV, 

M w = 0.78GeV, m v = 547.5MeV, m n > = 957.8MeV, 
M B = 5.2792GeV, r B o = 1.54 x l(T 12 s. (58) 

For the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, here we adopt the 
Wolfenstein parametrization for the CKM matrix up to O (A 5 ) with the parameters 
A = 0.2272, A = 0.818, p = 0.221 and fj = 0.3400. 

From Eq. fl45p we find the numerical values of the corresponding form factors at zero 
momentum transfer: AQ^ u (q 2 = 0) = 0.35 ± 0.05(co>fe) which are consistent with those 



given in 32 . 



B. Branching ratios 



For B° — > urf'' and 4>r)^ decays, the decay amplitudes as given in Eqs. (|54|) and (1561) 
can be rewritten as 



M 



V* b V ud T - V t * b V td P = V* b V ud T [1 + ze^} 



where 



vsy* 




p 


v: b v ud 




T 



a = art 



V u dV* b \ 



(59) 



(60) 



are the ratio of penguin to tree contributions and the weak phase (one of the three CKM 
angles) respectively, and S is the relative strong phase between tree (T) and penguin (P) 
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diagrams. The CP-averaged branching ratio for B° — > o;(0)r/^- ) decays can be then written 

as 

Br = (\M\ 2 + \M\ 2 )/2 = \V ub V: d T\ 2 [1 + 2z cos a cos 5 + z 2 ] , (61) 

where the ratio z and the strong phase 5 have been defined in Eqs.( l59l and ( |60i) . 

Using the wave functions and the input parameters as specified in previous sections, 
it is straightforward to calculate the branching ratios for the four considered decays. We 
find numerically that 

Br{ B° -> cur]) = [2.T^i{u b ) ± 0.7(a) ± 0.3(o a )] x 10~ 7 , (62) 
Br{ B° -> unf) = [0.75+°;^(^)^;^(a)+^ 9 7 (a 2 )] x 1(T 7 , (63) 

The main errors are induced by the uncertainties of ujb = 0.4 ± 0.04 GeV, a = 100° ± 20° 
and a 2 = 0.115 ± 0.115, respectively. 

For the B° — > <f>rft> decays, the decay amplitudes involve only the penguin (P) diagrams, 
the branching ratios are: 

Br { B° - 4>r)) = [6.3l^(o; 6 )l^(m s )l^(a 2 )] x 10^ 9 , (64) 
Br ( B° - H) = [7.3!!;^W^(m s )^(a 2 )] x 10" 9 , (65) 

The main errors are induced by the uncertainties of uj b = 0.4 ± 0.04 GeV, m s = 130 ± 30 
MeV and a 2 = 0.115 ± 0.115, respectively. 

For the CP-averaged branching ratios of the considered four decays, the PQCD predic- 
tions agree within errors with both the experimental upper limits as shown in Eqs.(JI]l2]), 
and the theoretical predictions in the QCDF approach, for example, as given in Ref. p|: 

Br( B° -> cor]) = (0.31+°,;^) x 10" 6 , (66) 
Br( B° -> ojt)') = (0.20lg;?J) x 10" 6 , (67) 
Br( B° -> 0r/ (/) ) w 1 x 10" 9 , (68) 

where the individual errors as given in Refs. [l| have been added in quadrature. It is 
easy to see that the central values of the pQCD predictions for Br[B — > (jjrft') are larger 
than the corresponding QCDF predictions, although they are still consistent if the large 
theoretical uncertainties are taken into account. The branching ratios Br{B —> ur]^) at 
the 10" 7 level can be measured in the forthcoming LHC experiments (33[. But it is very 
difficult to measure B — ► (fyqW decays because of their tiny (~ 10" 9 ) branching ratio, if 
the pQCD predictions are true. 

It is worth stressing that the theoretical predictions in the pQCD approach have large 
theoretical errors induced by the still large uncertainties of many input parameters. In 
our analysis, we considered the constraints on these parameters from analysis of other 
well measured decay channels. From numerical calculations, we get to know that the 
main errors come from the uncertainty of uu b , m s , a and a 2 ?w . Additionally, the final- 
state interactions remains unsettled in pQCD, which is non-perturbative in nature but 
not universal. 

In Fig. [3]we show the parameter dependence of the pQCD predictions for the branching 
ratios of B — ► lot] and iorj decays for u b = 0.4 ± 0.04 GeV, a = [0°, 180°]. From the 
numerical results and the figures we observe that the pQCD predictions are sensitive to 
the variations of u b and a. 
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FIG. 3: The a dependence of the branching ratios (in unit of 10 7 ) of B° — ► lot] and lot]' decays 
for (j) = 39.3°, u h = 0.36 GeV (dashed curve), 0.40 GeV (solid curve) and 0.44 GeV (dotted 
curve) . 

C. CP- violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of B — > corj^ decays 
in pQCD approach. Because these decays are neutral B meson decays, so we should 
consider the effects of B° — B° mixing. For B° meson decays into a CP eigenstate /, the 
time-dependent CP-violating asymmetry can be defined as 

Br f^(t)^f)-Br(B°(t)^f) 

= A d S r P cos(Amt)+A^sm(Amt), (69) 



Br (B\t)^f)+Br{B\t)^f) 



where Am is the mass difference between the two B® mass eigenstates, t — top — t tag 

is the time difference between the tagged B° (B ) and the accompanying B (B°) with 
opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct and 
mixing induced CP-violating asymmetries A^p and A^p can be written as 



(.;.,. \^cp\ 2 ~ 1 Amix _ 2Im(X C pi . 

Acp ~ i+\x CP r cp ~ i+\\ CP r { } 



12' ' I i \ ,.|2 

where the CP- violating parameter Xqp is 



_ VMf\Heff\^) _ « tt l + * e '<'-«> . . 

ACP V tb V t * d (f\H eff \B°) 1 + [n) 

Here the ratio z and the strong phase S have been defined previously. In pQCD approach, 
since both z and 5 are calculable, it is easy to find the numerical values of A% r P and A™™ 
for the considered decay processes. 
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a (degree) 



a (degree) 



FIG. 4: The direct and mixing-induced CP asymmetry (in percentage) of B° — ► ^(solid curve) 
and wf/ (dotted curve) decay as a function of CKM angle a. 

By using the central values of the input parameters, one found the pQCD predictions 
(in unit of 10~ 2 ) for the direct and mixing induced CP-violating asymmetries of the 
considered decays 

Af P {B* -> wn) = -69.li£i(a), A^(B° -> aj,) = +66.9^ 5 3 8 (a), (72) 
Af P {B« - = +13.9l«(a), - wi/) = +65.8^», (73) 

where the dominant errors come from the variations of a = 100° ± 20°. 

In Fig. HJ we show the a— dependence of the pQCD predictions for the direct and the 
mixing- induced CP- violating asymmetry for B° —>■ urj (solid curve) and B° — > urf (dotted 
curve) decay, respectively. The pQCD predictions in Eqs. ( 17211731) are consistent with the 
QCDF predictions due to very large uncertainties in the QCDF approach 

If we integrate the time variable t, we will get the total CP asymmetry for B° — > u)rft> 
decays (in units of 1CT 2 ): 

A&i&^ur,) = -ll±lti(a), (74) 
A%UB° - u»f) = +40.6™(a). (75) 

For the B° — > <fnf^ decay modes, however, there is no CP asymmetry, the reason is that 
they involve only penguin contributions and one type of CKM elements. 

It is worth of mentioning that although the CP-violating asymmetries of B — > lot)^ 
decays are rather large in size, but it may be difficult to measure them because of the 



predicted small branching ratios at 10 ' level 33 
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D. Effects of possible gluonic component of rf 



Up to now, we have not considered the possible contributions to the branching ratios 
and CP-violating asymmetries of B — > uj{<j))ri^ decays induced by the possible gluonic 
component of r\ and rf meson. For r/ meson, one generally believe that its gluonic com- 
ponent is negligibly small. For rf meson, however, a large gluonic component is generally 
expected because of the observed very large B — > Krf branching ratio. 

When a non-zero gluonic component exist in rf meson, an additional decay amplitude 
Ai' will be produced. Such decay amplitude may construct or destruct with the ones 
from the qq (q = u, d, s) components of rf , the branching ratios of the decays in question 
may be increased or decreased accordingly. 

In Ref. jof, by employing the pQCD factorization approach, Charng, Kurimoto and Li 
calculated the flavor-singlet contribution to the B — > 77W transition form factors induced 
by the Feynman diagrams with the two gluons emitted from the light quark of the B 
meson (see Fig. 1 of Ref. |9|]), the dominant gluonic contribution is negligible for B — > r\ 



form factor, and also small (less than 5%) for the B — > rf ones [13[ . 

Following the procedure of Ref. [9| and using the formulae given there, we also calculate 
the gluonic contributions to B — > uj(4>)r]' decays. If we add the gluonic contribution to 
the form factor F but keeping all other inputs in their default central values, we find 
numerically that 

Br( B° -> urf) = 0.82 x KT 7 , (76) 
Br( B° -> cprj') = 6.5 x KT 9 . (77) 

One can see that the variation of the central values of the pQCD predictions induced by the 
inclusion of gluonic contribution is only about ±10%, much smaller than the theoretical 
uncertainties from the variations of input parameters Ub, rn s , or a. 



V. SUMMARY 

In this paper, we calculated the branching ratios of B° — > ujr}^\ (jyq^ decays and CP- 
violating asymmetries of B° — > utjW decays at the leading order by using the pQCD 
factorization approach. 

Besides the usual factorizable diagrams, the non-factorizable and annihilation diagrams 
are also calculated analytically. Although the non-factorizable and annihilation contribu- 
tions are sub-leading for the branching ratios of the considered decays, but they are not 
negligible. Furthermore these diagrams provide the necessary strong phase required by a 
non-zero CP-violating asymmetry for the considered decays. 

After calculating all the diagrams, we found the branching ratios of the four related 
decays are very small, Br(B — > urf 1 ') are at the order of O(10~ 7 ) and Br(B — > <fnf'>) 
are around 10~ 9 . We also predict the direct and mixing-induced CP asymmetries of the 
B u)rf'' . Moreover, we compare our results with the current experimental data and the 
theoretical predictions in QCDF approach. In short, we found that 

• For the CP-averaged branching ratios of the considered decay modes, the pQCD 
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predictions are 



Br{ B° -> wrj) = (2.7l};J) x 10~ 7 , (78) 

Sr( 5° - wi/) = (0-75i°| 3 7 ) x lO" 7 , (79) 

Br( 5° -> 0r?) = (6.3±H) x 10~ 9 , (80) 

Br( 5° -> = (7.3±|:|) x 10~ 9 , (81) 

where the various errors as specified in Eqs. (|62 l) -( !65l) have been added in quadrature. 
The inclusion of the gluonic contribution can change the branching ratios of B — > 
lot}', <prf) decays by about 10% in magnitude. 

The pQCD predictions for the CP- violating asymmetries of B — > urj^ decays are 
large in size. 

The major theoretical errors are induced by the uncertainties of the input parame- 
ters u>b, the mass m s , the CKM angle a and the Gegenbauer moment a 1 ^ ■ 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function hi appearing in the expressions of the decay amplitudes in 
section Hill coming from the Fourier transformations of the function H^°\ 

h e (xi,x 3 ,b 1 ,b 3 ) = K (^/x^m B b 1 )[9(b 1 -b 3 )K (y/x^m B b 1 ) I (y/x^m B b 3 ) 

+9(b 3 - b\)K Q {y/x^m B b 3 ) I (y/x^m B bi)] S t (x 3 ), (Al) 

h a (x 2 , x 3 , b 2 , b 3 ) = K (iy/x 2 x 3 m B b 2 ) [9{b 3 - b 2 )K (iy/x^m B b 3 ) I (iy/x^m B b 2 ) 

+6{b 2 - b 3 )K {i^/x~ 3 m B b 2 ) I {i^ 3 m B b 3 )} S t (x 3 ), (A2) 

hf(x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = j#(& 2 - 6i)Io(M Bv /xix 3 6i)Ko(Mbv/^i^3&2) 

^ / K (M B F (1) b 2 ), for i? 1} >0\ 
+ (h ~ b 2 )j ■ ^ H W(M Bv %| b 2 ), for < oj ' ^ 

/ij(xi,x 2 ,x 3 ,6i,6 2 ) = - b^Kotiy/x^shM^Iotiy/x^^Ms) + (6 X «-> b 2 )\ 

■K (Vx[ + x 2 + x 3 - xix 3 - x 2 x 3 b x M B ), (A4) 
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h f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = < 9(h - b2)K (iy/x2X 3 b 1 M B )l (i^X2X 3 b 2 M B ) 



+ (h 



K^MsF^h), for Ff 2) > 



f kU(M b J\F*\ 60, for <0 



, (A5) 



(2) 



(2) 



where Jo is the Bessel function, K and Jo are the modified Bessel functions K (—ix) 
-(tt/2)Y (x) +i(ir/2)J (x), H^\z) = J (z) + iY (z), and the F {j) are defined by 



F (2) 



(x 1 - x 2 )x 3 
(xi - x 2 )x 3 



The threshold resummation form factor St(x-i) is adopted from Ref. [24 

2 1+2 T(3/2 + c 



v^r(i + c ) 



[x(l-x)] f 



where the parameter c = 0.3. This function is normalized to unity. 
The Sudakov factors appearing in section II III are defined as 



S ab (t) 



s cd (t) 



s (xim B / y/2, bA + s (x 3 m B / y/2, 63) + s ((1 - x 3 )m B / y/2, h 
1 



In ln(tM) , I In ln(t/A) 



-ln(6xA) -ln(M-) 
s {x\m B /y/2, b\ j + s [x 2 m B j y/2, b 

+s (x 3 m B / y/2, bA + s((l- x 3 )m B / y/2, b x 
1 



s (1 - x 2 )m B N2,b 



21nim + ln. ln ^ 



-ln(6iA) -ln(6 2 A)_ 
s ^iim B / y/2, 61 j + s (x 2 m B / y/2, b 2 ^j + s ^(1 - x 2 )m B / y/2, b 2 J 



+s (x 3 m B / y/2, 6 2 J + s ((1 - x 3 )m B / y/2, b 
1 



lnim + 21n. ln ^ 



-V(t) 



ln(6iA) -ln(6 2 A) 
s (x 2 m B / y/2, b 2 j + s (x 3 m B / a/2, 63 j + s f (1 - x 2 )m B / y/2, k 



+ 8 1(1 



x 3 )m B /y/2,b 3 ) - — 
/ pi 



In ln(t / A) , I In ln(t/A) 



- ln(6 2 A) 



ln(6 3 A)J 



(A6) 
(A7) 



(A* 



(A9) 



(A10) 



(All) 



(A12) 



where the function s(g, 6) are defined in the Appendix A of Ref. 19j. The scale t^s in the 
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above equations are chosen as 

t\ = max(y/x^m B , 1/h, l/b 3 ) 
t\ = max( A /x7m B , l/6 3 ) 
tl = max( v /x 3 "m jB , l/b 2 , 1/h) 
t\ = max( A /i^m B , l/b 2 , 1/b, 



3) , 



t f = max( v /a;iX3m i j, y/(xi - x 2 )x 3 m B , l/b u l/b 2 ) , 

t 3 f = max( v / xi + x 2 + x 3 - x\x 3 - x 2 x 3 m B , \Jx 2 x 3 m Bl l/&i, I/&2) , 

= max(v / xix 3 "m jB , - x 2 )x 3 m B , 1/bx, l/b 2 ) . (A13) 

They are given as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 
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